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Abstract 

We investigate the structure and stability properties of compact astrophysical objects that may 
be formed from the Bose-Einstein condensation of dark matter. Once the critical temperature 
of a boson gas is less than the critical temperature, a Bose-Einstein Condensation process can 
always take place during the cosmic history of the universe. Therefore we model the dark matter 
inside the star as a Bose-Einstein condensate, which can be described by a polytropic equation 
of state. We derive the basic general relativistic equations describing the equilibrium structure 
of the condensate dark matter star with spherically symmetric static geometry. The structure 
equations of the condensate dark matter stars are studied numerically. The critical mass and 
radius of the dark matter star are given by M cr a ~ 2(7 a /l/m.) 1 / 2 (m x /l GeV)~ 3 / 2 AfQ and R cr it ~ 
1.1 x 10 6 (7 a /l fm) 1 / 2 (m x /l GeV)~ 3 / 2 cm respectively, where l a and m x are the scattering length 
and the mass of dark matter particle, respectively. 
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I. INTRODUCTION 



Since in the terrestrial experiments Bose-Einstein Condensation (BEC) is a well-known 
phenomenon, the possibility that a similar condensation may have occurred during the 
cosmo logical evolution of the universe cannot be excluded a priori. In fact, once the critical 
temperature of the boson gas is less than the critical temperature, BEC can always take 
place at some moment during the cosmic history of the universe. Different aspects of the 
BEC cosmological transition were considered in The critical temperature for the 

condensation to take place is T cr < 27rh 2 n 2 ^ 3 /m x kB, where n is the particle number density, 
m x is the particle mass, and ks is Boltzmann's constant fl-js]. Since the matter temperature 
T m varies as T rn oc a -2 , where a is the scale factor of the universe, it follows that during an 
adiabatic evolution the ratio of the photon temperature T 7 and of the matter temperature 
evolves as T 7 /T m oc a. Cosmic evolution has the same temperature dependence, since in an 
adiabatic expansion process the density of a matter dominated universe evolves as p oc T 3 / 2 
{2I, 3|. Therefore, if the boson temperature is equal, for example, to the radiation temperature 
at a redshift z = 1000, the critical temperature for the Bose-Einstein Condensation is at 
present T cr = 0.0027-K" [3]. On the other hand, we expect that the universe is always under 
critical temperature, if it is at the present time j^]. 

It has been proposed recently that the first stars to exist in the universe were powered 
by dark matter heating rather than by fusion j^-Q. Weakly Interacting Massive Particles, 
collect inside the first stars and annihilate, to produce a heat source that can power the 
stars. A new stellar phase results, called a Dark Star, powered by dark matter annihilation 
as long as there is dark matter fuel. The heat source can power the Dark Stars for millions 
to billions of years. These objects can grow to be supermassive dark stars (SMDS) with 
masses ~ (10 5 — 1O 7 )M J9]. However it is not known what is the ratio between dark 
matter and anti-dark matter or if the dark matter particle should even have an anti-particle 
partner to annihilate to generate the heat. On the other hand non- annihilating dark matter 
particles can also generate very important effects in main sequence stars. In lOj it was 
shown that the energy transport mechanism induced by dark matter particles can produce 
unusual conditions in the core of main sequence stars, which constrain the spin-dependent 
cross section > 10 _37 cm 2 , and the dark matter particle mass m x > 5GeV. 

It is the purpose of the present paper to consider a systematic study of the effects of the 
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presence of condensate dark matter on the properties of compact astrophysical objects. By 
considering the condensation process in a cosmological setting, as a first step in our analysis 
we determine the physical and cosmological parameters at which the condensation of a boson 
gas can take place. If the critical temperature of the gas is less than the critical temperature, 
the Bose-Einstein Condensation process takes place, resulting in the formation of a conden- 
sate dark matter background, which can be either accreted by normal matter baryonic stars, 
or can form stellar type objects due to gravitational instabilities. Therefore we model the 
dark matter inside the star as a Bose-Einstein condensate. The condensate dark matter 
equation of state can be described by a polytropic equation of state, with polytropic index 
n = 1. We derive the basic general relativistic equations describing the equilibrium structure 
of the condensate dark matter star with spherically symmetric static geometry, and we study 
the structure equations numerically. Our results show that the critical mass and radius of 
the condensate dark matters star are given by M crit « 2(/ a /l/m) 1 / 2 (m x /l GeV)~ 3 / 2 M 
and Rent ~ 1.1 x 10 6 (Z a /l fm) 1//2 (m x /l GeV) -3 / 2 cm, respectively, where l a and m x are the 
scattering length and the mass of the dark matter particle. 

The present paper is organized as follows. In Section [III the equations of state of the 
condensate and non-condensate dark matter are written down. The structure equations 
describing the properties of the general relativistic dark stars are derived in Section IIII1 The 
results of the numerical integration of the structure equations are presented in Section IIV1 
We discuss and conclude our results in Section [V] 



II. DARK MATTER EQUATIONS OF STATE 
A. Bose-Einstein condensed dark matter 



At very low temperatures, all particles in a dilute Bose gas condense to the same quantum 
ground state, forming a Bose-Einstein Condensate (BEC). Particles become correlated with 
each other when their wavelengths overlap, that is, the thermal wavelength is greater than 
the mean inter-particles distance /. This happens at a temperature T cr 2%xh 2 p 2 / 3 /m 5 / 3 kB, 
where m is the mass of the particle in the condensate, p is the density, and is Boltzmann's 
constant 4j. A coherent state develops when the particle density is enough high, or the tem- 
perature is sufficiently low. We assume that the dark matter halos are composed of a strongly 



3 



- coupled dilute Bose-Einstein condensate at absolute zero. Hence almost all the dark mat- 
ter particles are in the condensate. In a dilute and cold gas, only binary collisions at low 
energy are relevant, and these collisions are characterized by a single parameter, the s-wave 
scattering length l a , independently of the details of the two-body potential. Therefore, one 
can replace the interaction potential with an effective interaction Vj (r ' — f) = X5 (r ' — f), 
where the coupling constant A is related to the scattering length l a through A = A7rh 2 l a /m x 
4j. The ground state properties of the dark matter are described by the mean-field Gross- 
Pitaevskii (GP) equation. The GP equation for the dark matter halos can be derived from 
the GP energy functional, 

h 2 
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where ip (r) is the wave function of the condensate, and Uq = Anh 2 l a /m x \&\. The first term 
in the energy functional is the kinetic energy, the second is the interaction energy, and the 
third is the gravitational potential energy. The mass density of the condensate dark matter 
is defined as 

Px ( f ) = m x 1^ (Ol 2 = m xP *) > ( 2 ) 
and the normalization condition is N = J \ip (r)\ 2 df, where is the total number of dark 



matter particles. The variational procedure SE 
equation as 

h 2 
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V ip (r) + m x V (f) tjj (r) + Uq (r) | ip{f) — flip (f) 



(3) 



where /i is the chemical potential, and the gravitational potential V satisfies the Poisson 
equation 

V 2 V = AirGp. (4) 

In the time-dependent case the generalized Gross-Pitaevskii equation describing a gravita- 
tionally trapped rotating Bose-Einstein condensate is given by 
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(5) 



The physical properties of a Bose-Einstein condensate described by the generalized Gross- 
Pitaevskii equation given by Eq. (jSJ) can be understood much easily by using the so-called 



Madelung representation of the wave function |4|-|6|], which consist in writing ip in the form 



(r,t) = a/ p (f, t) exp 



h 



S(r,t) 



(6) 



where the function S (r, t) has the dimensions of an action. By substituting the above 
expression of i/j (f, t) into Eq. (JSJ), it decouples into a system of two differential equations for 
the real functions p x and v, given by 



dt 



+ V • ( Px v) = 0, 



(7) 
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where we have introduced the quantum potential Vq — — (h 2 /2m x ) V 2 Ai /p^/ v /p^, and the 
velocity of the quantum fluid v = VS/m x , respectively The effective pressure of the con- 
densate is given by 

(9) 
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The Bose-Einstein gravitational condensation can be described as a gas whose density 
and pressure are related by a polytropic equation of state, with index n = 1 {4, fj. When the 
number of particles in the gravitationally bounded Bose-Einstein condensate states becomes 
large enough, the quantum pressure term makes a significant contribution only near the 
boundary of the condensation. Thus the quantum stress term in the equation of motion 
of the condensate can be neglected. This is the Thomas-Fermi approximation, which has 
been extensively used for the study of the Bose-Einstein condensates 4|t6[- As the number 
of particles in the condensate becomes infinite, the Thomas-Fermi approximation becomes 
exact. This approximation also corresponds to the classical limit of the theory. From its 
definition it follows that the velocity field is irrotational, satisfying the condition V x v — 0. 



B. Normal dark matter equation of state 

We assume that in the early stages of the evolution of the universe dark matter consisted 
of bosonic particles of mass m x and temperature T, originating in equilibrium and decoupling 
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at a temperature To or chemical potential fi » m x . By assuming that the dark matter 
forms an isotropic gas of particles in kinetic equilibrium, the spatial number density is given 
by 

where h is Planck's constant, g is the number of helicity states, and 

f(p) = [ex V (E-fi)-ir\ (12) 

where p is the momentum of the particle and E = ^Jp 1 + m^c 4 is the energy. A particle 
species that decouples in the early universe from the remaining plasma at temperature Tp 
redshifts its momenta according to p(t) = Pd^d/ 'a{t), where a(t) is the cosmological scale 
factor and od is the value of the scale factor at the decoupling. The number density of 
the particles n evolves as n x ~ a~ 3 (t) 12J. The distribution function / at a time t after 
the decoupling is related to the value of the distribution function at the decoupling by 
f(p) = f (pa/a,E)). The distribution function keeps an equilibrium shape in two regimes. In 
the extreme- relativistic case, when E ~ pc, T = Ti>ao/ci, and fi = findo/a, respectively, 
the distribution function is given by /er(p) = [exp (pc — jj) — 1\~ . In the non- relativistic 
decoupling case E — \i « p 2 /2m x — fi^ n , where we have defined /ikin = A* — m x c2 ' ^he 
distribution function is /nr(p) = [exp(p 2 /2m x — jj kill ) — 1] . In the non-relativistic case 
/ikm and T evolve as //kin = A^kin.D i a D/ a ) 2 and T = T D / (a D /a) 2 , respectively 12]. 

The kinetic energy-momentum tensor associated to the frozen distribution of dark 
matter is given by 



T '-4?l^r ! ^ < 13) 



The energy density e of the system is defined as 

e = wJ Ef ^ d ^ ^ 

while the pressure of a system with an isotropic distribution of momenta is given by 

p = ~h / pvf ^ d3p = wj irrtp)^' ( 15 ) 

where the velocity v is related to the momentum by v = pc 2 /E [13] . In the non-relativistic 
regime, when E « m x c 2 and p ~ m x v x , the density p x of the dark matter is given by 
Px = m x n xi wn ile its pressure P x can be obtained as [12] 



m x 
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giving 

Px = Px<?°l (17) 

where a 2 = {v 2 )/3c 2 , and (v 2 ) is the average squared velocity of the particle. a v is the 
one-dimensional velocity dispersion. 

C. The cosmological Bose-Einstein transition process 

In order to analyze the cosmological conditions for the formation of a dark matter Bose- 



Einstein condensate we follow the approach introduced in [14j . The cosmological parameters 
at which the Bose- Einstein condensation process took place can be estimated by taking into 
account that the laws of thermodynamics require that both the chemical potential p and the 
pressure p are single valued functions, that is, for any given values of the particle number n 



and T, there must only exist a single value of p or p, respectively [15]. 

Therefore, a first thermodynamic condition that must be satisfied during the cosmological 
Bose-Einstein Condensation process is the continuity of the pressure at the transition point. 
With the use of Eqs. (Q and (IT7|) the continuity of the pressure uniquely fixes the critical 
transition density p° r from the normal dark matter state to the Bose-Einstein condensed 
state as 

The numerical value of the transition density depends on three unknown parameters, 
the dark matter particle mass, the scattering length, and the dark matter particles velocity 
dispersion, respectively. By assuming a typical mass of the dark matter particle of the order 
of 1 GeV (1 GeV = 1.78 x 10~ 24 g), a typical scattering length of the order of 1 fm, and a 
mean velocity square of the order of (v 2 ) = 9 x 10 14 cm 2 /s 2 , the critical transition density 
can be written as 

rr 2 \ / m \ 3 / 7 \ _1 



by 



£ = 7.327 x 10" [j^)^) {^) g/cm<. (19) 
?he critical temperature at the moment of Bose-Einstein condensate transition is given 
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7 



where C(3/2) is the Riemann zeta function, or 

r ^ 3 - 76xl08 (r^v) 1/3 (i^) !/3 (ii)" 2/3K - < 21 » 

The critical pressure of the dark matter fluid at the condensation moment can be obtained 

as 

- "i x ^ (l^) 2 (l^v) 3 (t^)" 1 < 22 > 

By assuming that the Universe is flat, the cosmological dynamics of the dark matter 
before the condensation is described by the standard Friedmann-Robertson- Walker model, 
in which the evolution of the normal dark matter density is given by the equation 

p x + 3p x (l + a 2 )- = 0, (23) 

(X 

where a is the scale factor of the Universe, with the general solution given by 



where p Xi o is the density of the dark matter at a = a . Usually the scale factor is normalized 
so that at the present time a — ao — 1. Therefore the critical value a cr of the scale factor 
of the Universe at the moment of the beginning of the Bose-Einstein condensation can be 
immediately obtained as 



<-/<» = = w g l (25) 



where we have introduced the critical density of the Universe p cr> o = 3Hq/StiG, and the dark 
matter density parameter Q x0 = p Xi o/p C r,o ? where H is the present day Hubble constant. 
Hence for the value of the cosmological redshift at which the Bose-Einstein transition did 
occur we obtain the expression 



J-xfi 

In the following for the Hubble constant we adopt the value i^o = 70 km/s/Mpc 



2.27 x 10 18 s l , giving for the critical density a value of p crj o = 9.24 x 
present day dark matter density parameter is f2 x ~ 0.228, respectively 



q-30 g / cm 3_ The 



16| . By using the 
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adopted numerical values of the constants we obtain for the critical scale factor and for the 
critical redshift at which the Bose-Einstein condensation took place the values 



a 

and 



. / a 2 \ -1/3(1+^) , x Nl/3(l +CT 2 ) 



/m \(i+- 2 )/ a 2 \ 1/3 ( 1+<j2 ) / / x-i/3(i+- 2 ) 
respectively. 

These results show that if dark matter consists of self-interacting particles with masses 
of the order of 1 GeV, and scattering length of the order of 1 fm, the transition from normal 
dark matter to condensate dark matter did occur very early in the history of the Universe, 
presumably even during the post-inflationary reheating phase. These results are consistent 
with those obtained in |3| , where the inflation is naturally initiated by the condensation of the 
bosons in the huge vacuum energy. Bose-Einstein condensation can take place and continue 
provided the boson temperature is less than the critical temperature at some moment of 
cosmic evolution. The corresponding critical boson mass for BEC can be arbitrary in general, 
and there are no theoretical restrictions for its value jsj]. After the beginning of the phase 
transition the density of the dark matter p x (t) decreases from (T cr ) = p£ or (when all the 
dark matter is in a normal, non-condensed form) to p x (T cr ) = p x EC , corresponding to the 
full conversion of dark matter into a condensed state. 



III. HYDROSTATIC EQUILIBRIUM EQUATIONS OF DARK MATTER STARS 

In the following we restrict our study of the condensate dark matter star to the static 
and spherically symmetric case, with the metric represented as 

ds 2 = e^c 2 dt 2 - e x ^dr 2 - r 2 {dd 2 + sin 2 M<f) 2 ) . (29) 

For the metric given by Eq. ( |29|) . the Einstein gravitational field equations, describing 
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the dark matter star, take the form 



A'\ 1 8nG , . 



,V 1\ 1 8ttG 
e" A T + 15 " 15 = ( 31 ) 



2 / ,y 2 
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2 _ V 2 r 2 J c 4 

and 



0"' 



i/ = -2 , (33) 

where ' = d/dr, e is the total energy-density, and a is the pressure along the radial direc- 
tion. Equation (1331) is the consequence of the conservation of the energy-momentum tensor, 
T^u-^ = 0. Eq. (1301 can be easily integrated to give 

\ 2GM(r) , . 

e- A = 1 ^1, 34 

where M(r) = 4n J (e/c 2 ) r 2 dr. With the use of Eqs. (I3TI) and (|33|) we obtain the mass 
continuity equation, and the generalized Tolman-Oppenheimer-Volkoff (TOV) equation de- 
scribing the spherical symmetric static dark matter configurations, 

— = 4vr-r 2 , 35 
dr cr 

da (e + a) [(47rG/c 4 ) ar 3 + GM/c 2 ] 

d^~ r 2 (1 - 2GM/c 2 r) ' ^ ' 

In the case of the dark matter stars e = p x c 2 and a = P x = (2iih 2 l a /'m 3 x j p 2 ,. In order 

to simplify the equations we introduce the dimensionless dark matter density 8 X , defined as 

Px = Pcx®x-> where p cx is the central density of the dark matter. 

In addition, let us denote 

uppcx 2nh 2 l a / l a \ ( m x \~ 3 ( p cx \ 

a x = —r- = -^Pcx = 1-658 — - — in1 . , —z , (37) 



c 2 c 2 m 3 \lfmy \m p J \10 16 g/cirr 

The coefficient a x can be broken into two parts. Let a x = /3p cx , where 



2Txh 2 l a . _l« / / 



= ^fi- = 1.658 x 10" 16 | ~r~r ) ( — — J cm 3 /g. (38) 
c z m x \ 1 1m / \ m p J 

(3 is only determined by the properties of dark matter particle. Furthermore, a x can be 

constrained by the condition that the sound speed of condensate dark matter must be 

smaller than the speed of light. Since 

P x = ^c 2 pl (39) 

Pc X 
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the constraint on sound speed < c gives a x < This condition must be valid 

everywhere inside the star, at the center of the center it gives 

a x < 1 -. (40) 

Later we will find that a xmax = 0.43 from the instability analysis. 
By introducing a set of dimensionless variables (rj, m), defined as 

V, (41) 



and 

1 /c 2 \ 3/2 1 



M — —= — -— m, (42) 



the structure equations of the dark matter star can be written as 

x x dr] r) 2 (l- 2m/r]) ' 1 ' 

The initial conditions for the structure equations are 

m(0) = 0, (45) 

X (0) = 1. (46) 

And boundary of the object is reached when 9 X = at which point the values for m and r\ 
are denoted by m and 770- 



IV. NUMERICAL RESULTS 



The solutions of the dimensionless equations (143 j) and (I44p are studied numerically for 
some fiducial values m x = m p , l a = 1 fm, and p cx = 10 13 g/cm 3 . The variation of the 
dimensionless density 9 X as a function of the dimensionless radius 77 is plotted in Fig. [TJ 

The only parameter in the structure equations of the star Eqs. (143]) and ( )44|) is a x . The 
dimensionless mass m and the dimensionless radius r are thus functions of a x only. The 
dependence of m and 770 on a x can be fitted by simple relations valid for physical choice of 
a x < 1/2 

m = 8.89a 3 / 2 (1 + 6a x )~ 4/3 , (47) 
11 
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FIG. 1: Dimensionless density 9 X versus dimensionless radius rj for condensate dark matter stars 
with different values of a x . 

7/o = 4.45aJ /2 (1 + 6a x )" 2/5 . (48) 

Figs. |2] and [3] illustrate the comparison between the numerically obtained values and the 
fitted values as functions of a x . Fig. H] illustrates the relative error of the fitted values from 
Eqs. (|4T|) and (148!) . and it shows that they are very good approximation formulae. It can be 
seen that Eq. ( ]4"T|) and Fig. HE] give very good fitting to the numerically obtained values. 

Numerical m 

10° -a Fitted m 




10^ 10" 5 10-' 10 J W 2 10 - ' 



FIG. 2: Comparison between the numerically obtained mo and the fitted mo from Eq. (|47|) : the 
black line is for the numerically obtained mo, and the red line gives the fitted mo from Eq. (|47p . 
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FIG. 3: Comparison between the numerically obtained rjo and the fitted r/o from Eq. f|48|) : the 
black line represents the numerically obtained rjo, and the red line represents the fitted i]q from 
Eq. (0H]). 




FIG. 4: Relative error of the fitted mo and 770 from Eqs. (|47p and (|48p : the black line represents 
mo and the red line represents 770 . 

The dimensionless mass mo increases with the dimensionless radius rjo, and there is no 
instability. However, when Eq. ( 142 jl and Eq. ( HIT) are used to calculate the physical mass and 
radius, they will behave quite differently. When the central density keeps increasing, there 
is a maximum mass corresponding to the instability. The instability is clearly seen in Fig. [5] 
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for different choices of /3 when dM/dR = 0. We will explain the origin of the instability in 
the next Section. Fig. [6] illustrates the dependence of the physical mass on the parameter 
a x = (3p cx for different values of j3. 

There is a simple criteria for instability onset 



Pp, 



0.43, 



(49) 



which is smaller than the constraint value obtained by the sound speed limit. In the next 
Section we will explain the origin of this value. 

By using Eq. (l47j) and Eq. (|48l) . we can derive simple relations between the critical mass 
and radius corresponding to the instability, and the dark matter particles property parameter 



P, 



and 



M Q ( —J- 3 
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(50) 



R crit = 8.75 x 10 5 cm 



P 



10- 16 cm 3 /g 



1.1 x 10 6 (l a /l fm) 1/2 (m x /l GeV)- 3/2 cm, (51) 



respectively. 
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FIG. 5: Numerical relation between physical mass and radius for different choices of /3: the dashed 
lines correspond to non-physical condition where a x > 1/2. 
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FIG. 6: Numerical relation between physical mass and a x for different choices of j3: the dashed 
lines correspond to non-physical condition where a x > 1/2. 

V. DISCUSSIONS AND FINAL REMARKS 

In the present study, we have shown that it is possible to form stable dark matter stellar 
objects, given that dark matter has self-interaction determining the condensation process. 
The mass and the radius of the star depend on two parameters, one being an "environmental" 
parameter, the central density p cx , and the other is the intrinsic property of dark matter 
particle (3. 

The parameter (3 is a combination of dark matter particle mass m x and the self-interaction 
scattering length l a , which is related to the self-interaction cross section a xx by the relation 
a xx = 47r/^. Currently, the most popular dark matter model, the Weakly- Interacted-Massive- 
Particle (WIMPs) model favors a mass range between GeV to several TeV. A possible dark 
matter candidate particles, the sterile neutrinos, have a mass in the keV range. This mass 
range cannot explain the cold dark matter structure formation. Other models of dark matter, 



3, 



ike supersymmetric particles, light gravitons, axions predict different mass ranges (cf. 

1^| )• The idea that dark matter may have self interaction was first proposed by Spergel and 
Steinhardt 19j to alleviate several apparent conflicts between astrophysical observations 
and the collisionless CDM model. Their proposed self-interaction has the strength range 
<T xx/ rn x = 0-45 — 450 cm 2 /g. The experimental constraint of self-interaction cross section 
mostly come from astronomical observation of colliding Bullet Cluster with the upper limit 
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<T xx/ m x < 1 cm 2 /g 20j. The condensate dark matter star might be able to provide another 
means to constrain the properties of dark matter particles. 

The simple criteria for instability, Eq. ( I49p . might have a simple explanation by calculating 
the total energy of the condensate dark matter star. The total energy is the sum of negative 
gravitational energy E gTSLV and the internal energy E [nt . By Newtonian approximation the 



gravitational energy is 



g g ,v - I -g^dm = - G jl Gpex f«*xVbl, (52) 
and the internal energy of the condensate dark matter is 

E -=S u < AV=a *^mzl e2 ^ (53) 

At the instability point, E = E gTaiV + E int = 0. This simplifies to be 

a x J O^fdrj — J m6 x r]dr] = 0. (54) 

Since the two integrals are only functions of a x , the instability onset corresponds to the zero 
point of an equation of a x , which leads to Eq. (|4"9"1) . 

When a condensate dark matter is formed, it will lead to the accretion of material from 
space. The incoming dark matter particles will interact with the dark matter in the con- 
densate star, transfer their kinetic energy, and get captured by the condensate star. The 
accretion process will increase the mass of the condensate star, and may exceed the instabil- 
ity limit, thus causing the collapse of the condensate star. This collapse will probably form 
black holes from the dark matter star. 

In addition to the formation of the dark matter stars during the phase transition epoch 
in the early universe, dark matter stars may also be formed in another ways. For example, 
the formation of early stars will make a deep gravitational well to capture the dark matter, 
which has a high density in the early universe. This process may result in a condensate 
dark matter core to form inside normal stars. As more dark matter is captured, due to the 
self-interaction between particle the accretion rate will be enhanced. The condensate dark 
matter core will collapse to form black holes once its mass exceeds the instability limit and 
swallow the degenerate core of normal matter core, if it exists. This collapse process might 
have signatures in the observations of high redshift long Gamma-ray Bursts. 
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As it was pointed out in |z-9], the annihilation of the dark matter particles inside the 
star can dramatically alter the evolution of the system, by providing a significant source of 
heating. The dark matter particles annihilation rate is n x < av >, where the standard anni- 
hilation cross section has values of the order of < av >~ 3 x 10 -26 cm 3 /s. The dark matter 
particle annihilation produces energy at a rate per unit volume of Qdm =< & v > P x / m x- 
The annihilation products typically are electrons, photons, and neutrinos. The neutrinos 
escape the star, while the other annihilation products are trapped in the dark star, and 
thermalize within the star, eventually heating it up. The luminosity from the dark matter 
particle annihilation is I/dm = Jq j Qdm<IV , where fq is the fraction of the annihilation 
energy deposited in the star. Due to their high density, the dark matter condensate stars 
have a high neutrino luminosity, and this intense neutrino flux may play a significant role 
in the formation of the gamma ray bursts. 

It is also interesting to compare condensate dark matter star to the non-condensate case. 
In Section IIC, the equation of state for non-condensate dark matter can be described by 
the equation of state of a collisionless system, P x = p x c 2 a%, where a v is the dimensionless 
velocity dispersion. This equation of state corresponds to the polytropic model with n = oo. 
Since the dark matter density in non-condensate state is much smaller than that of the 
condensate, the radius of the non-condensate dark matter star must be much larger than 
that of the condensate star. Therefore we can use the Lane-Emden equation to describe the 
stellar structure for the non-condensate case, i.e. 

= ,55) 

r z ar \p x ar J c z a z 

Eq. ( )55l) . describing the non-condensate case, corresponds to the isothermal sphere, and has 
an analytical solution given by 

ftfr) - (») 

This solution is unstable and leads to singular behavior of the density near the center. It 
also does not have a definite boundary with finite mass. This will lead to an extended 
object made of normal dark matter, while in the condensate case, the dark matter star is 
very compact. 

Alternatively, when normal dark matter is able to form stellar object, it might be de- 
scribed as an adiabatic process, P X P X ^ = const. In the simplest case, dark matter only has 
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3 degrees of freedom, i.e. 7 = 5/3, and the equation of state is given by 

P x = Kpf. (57) 
The constant K can be determined from the adiabatic invariants, 



-5/3 _ p -5/3 _ Pcx u rp -5/3 _ ksTcr 2/3 



"X ""X 

where P c and p cx are central pressure and the density, respectively. The central temperature 
is taken to be the phase transition temperature. The structure is thus described by the 
Lane-Emden equation for n = 3/2. The mass and radius are given by 

3/2 



M = 2.714 x 4tt (59) 



and 

R = 3-654 x^J p- 1 /-, (60) 

respectively. For the fiducial values m x = m p and l a — 1 fm, we obtain T cr = 3.68 x 10 8 K. 
The mass and radius are given by 

1/2 

lg/cm 3 y 

\ -1/2 

Pcx 



M = 476.15M ( Pcx ) , (61) 



i? = 15.82i? -924 . (62) 
\lg/cm d / 

From these equations it can be clearly seen that for a given mass non-condensate dark 
matter stars have a much higher radius. This result also justifies the use of the Newtonian 
Lane-Emden equation to describe this situation. 

If the Bose-Einstein condensation took place in the very early universe, at very high 
redshifts, as suggested by our analysis of the phase transition, at the moment of star/galaxy 
formation most, if not all, of the dark matter was already in a condensate form. Hence 
one may naturally expect the presence of the accreted condensate dark matter inside stellar 
type objects in the early stages of star and galaxy formation, as well as the formation, 
through gravitational instabilities, of pure condensate dark matter stars. Thus Bose-Einstein 
condensate dark matter stars may have been more common in the early Universe than normal 
dark matter stars. 
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